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Abstract 

We study the local approximation properties in hierarchical spline spaces through 
multiscale quasi-interpolation operators. This construction suggests the analysis of 
a subspace of the classical hierarchical spline space [VGJSll] which still satisfies 
the essential properties of the full space. The B-spline basis of such a subspace can 
be constructed using parent-children relations only, making it well adapted to local 
refinement algorithms. 
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1 Introduction 

Local adaptivity in numerical methods for partial differential equations makes possi¬ 
ble to solve real problems leading to a suitable approximation of the desired solution 
without exceeding the limits of available software. When considering isogemetric meth¬ 
ods [HCB05, CHB09], from a theoretical point of view, the design of efficient and robust 
strategies for local rehnement constitutes a challenging problem because the tensor prod¬ 
uct structure of B-splines [dVOl, S07] is broken. 

Hierarchical B-splines (HB-splines) based on the construction presented in [K98, K97, 
VGJSll] are a promising approach, because their construction and properties are closely 
related to the ones of hierarchical hnite elements. Moreover, truncated hierarchical B- 
splines (THB-splines) have been introduced in [GJS12], where their use as a framework 
for isogeometric analysis that provides local rehnement possibilities has been analysed; 
see also [GJS14]. Local approximation estimates for hierarchical spline spaces have been 
studied in [SM14] using quasi-interpolants described in terms of the truncated hierarchical 
basis while the use of THB-splines in conjunction with residual based error indicators has 
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been proposed in [BuGilS] under a few assumptions on the meshes. It is important to 
remark that truncation is indeed a possible strategy to recover partition of unity and 
convex hull property. On the other hand, the procedure of truncation requires a specihc 
construction that entails complicated basis function supports (that may be non convex 
and/or not connected) and their use may produce a non negligible overhead with an 
adaptive strategy. 

In the present paper, we take the point of view of classical hierarchical B-splines and 
we study their structure with a special attention to all those properties that may be 
needed or may facilitate their use with an adaptive isogeometric method. 

When considering an underlying sequence of nested tensor-product spline spaces and 
the corresponding B-spline bases, a particular way of selecting B-spline basis functions 
from each different level in order to build a hierarchical basis "H has been established 
in [K98] (see also [VGJSll]). This hierarchical basis l-L enjoys some important properties: 

• H is a set of linearly independent B-spline functions. 

• It is possible to identify uniquely the basis l-L from the knowledge of a hierarchical 
mesh. 

• All functions in the coarsest underlying tensor-product spline space belong to the 
hierarchical space spanH. 

• Under certain assumptions about the hierarchy of subdomains associated to H, it 
is possible to dehne a multiscale quasi-interpolant operator in span?/ obtaining 
optimal orders of local approximation. 

• Any enlargement of the hierarchy of the subdomains associated to H gives rise to 
a refined basis H* in the sense that spanH C spanH*. 

On the other hand, unlike tensor-product B-spline bases, the functions in the hierar¬ 
chical basis H do not constitute a partition of unity. If we consider the corresponding 
coefficients C M in order to form such a partition, i.e., 

^0/3/3 = !, (1) 

it is known that the coefficients are nonnegative, but in fact, some of them can be 
equal to zero. 

In this article we analyse the local approximation properties of hierarchical splines 
spaces through the construction of a multiscale quasi-interpolant operator. Kraft [K98] 
has introduced such a operator for the case of bivariate spline spaces on inhnite uni¬ 
form knot vectors and has studied its pointwise approximation properties. We extend 
his results to the case of open knot vectors with possible multiple internal knots in d- 
dimensional domains, for d > 1, and we also provide local approximation estimates in 
2,9-norms, for 1 < g < cx). 

Furthermore, we propose a new hierarchical spline space, through a construction of a 
set of basis functions named H that satishes all important properties just listed above. 
In particular, H G H, and therefore, in general, the new hierarchical space span?/ may 
be smaller than spanTf. However, the new basis is easier to manage and to update when 
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performing local adaptive refinement, and in this case, the coefficients for (1) are strictly 
positive. 

This paper is organized as follows. In Section 2 we introduce the notation and the 
assumptions for the underlying tensor-product spline spaces to be considered, and in 
Section 3 we briefly introduce the standard hierarchical B-spline basis and prove some 
results which will be useful later. In Section 4 we construct a multiscale quasi-interpolant 
operator and study the local approximation properties in hierarchical spline spaces. In 
Section 5 we define the new hierarchical B-spline basis and prove some of its basic prop¬ 
erties. Finally, we conclude the article with some final remarks in Section 6. 

2 Spline spaces and B-spline bases 

Univariate B-spline bases Let := be a p-open knot vector, i.e., a 

sequence such that 

0 ^ ^p+‘2 — ^ C^n+1 ■ ■ ■ ^n+p+l 1; 

where the two positive integer p and n denote a given polynomial degree, and the cor¬ 
responding number of B-splines defined over the subdivision respectively. Here, 
n > p + 1. We also introduce the set '■= {Cj}^=i of knots without repetitions, and 
denote by ruj the multiplicity of the breakpoint Cp such that 

'— 'p^n { Cl? • ; C2^ • • • ‘ ^ Cll {? 

mi times m2 times mu times 


with mi = n + p + 1. Note that the two extreme knots are repeated p -|- 1 times, i.e., 

i=l 

mi = ma = p -|- 1. We assume that an internal knot can be repeated at most p -|- 1 times, 
i.e., mj < p + 1, for j = 2,. .., fi — 1. 

Let B(Ep^n) ■= {^ 1 ? b 2 ,..., bn} be the B-spline basis (cf. [dVOl, S07]) associated to the 
knot vector The local knot vector of bj is given by 

^bj •= {Cj? • ■ ■ ? Ci+p+l}? 

which is a subsequence of p -|- 2 consecutive knots of Sp^„. We remark that 

supp6j = [CpCj+p+i]- 
Let X(Sp „) be the mesh defined by 

X(Sp,„) := {[CpCi+i] |j = 

For each I = [Cj, Cj+i] ^ ^{'^p,n) there exists a unique k = J2i=i ®^oh that / = [Cfc, Cfc+i] 
and Cfc 7 ^ Cfc+i- The union of the supports of the B-splines acting on / identifies the support 
extension /, namely 

• [Cfc—p? Cfc+p+i] • 
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Let So := be a given p-open knot vector. We consider a seqnence of 

sucessive refinements of Sq, i.e., 

SoCSiC..., (2) 

where is a p-open knot vector, and C means that Ei is a snbsequence of S^+i, 
for £ G Nq. In other words, condition (2) says that if ^ is a knot in with multiplicity 
m, then ^ is also a knot in S^+i with multiplicity at least m, for £ G Nq- 
Let Bi := B{Ei) be the B-spline basis, for i G Nq. 

Definition 2.1. Let Pi G Bi be given and let be the corresponding local knot vector. 
Let C EiJ^i be the knot vector obtained from E^^ after inserting the knots of 

which are in the interior of supp Pi. We say that G is a child of Pi G Bi if the 
local knot vector of Pi+i, which is denoted by is a subsequence of In other 

words, the children of Pi are the B-splines in whose local knot vector consist of p + 2 
consecutive knot of We let 

C{Pi) := {Pi+i G Bi+i I Pi+i is a child oi Pi]. 

Conversely, if G BiJ^i is given, we dehne the set of parents of Pi+i by 

V{Pi+i) := {Pi G Bi I Pi+i is a child of Pi}. 

It is easy to check that the last dehnition means that /S^+i G is a child of Pi G Bi 
if and only if 

(i) minS^^ < minS^^^^ < maxS;3^_^^ < maxS^^. 

(ii) If ^ G matches any of the end points of then the multiplicity of ^ in 
is less or equal to the multiplicity of ^ in Ej^^. 

In particular, notice that if /d^+i is a child of /?£, then supp /d^+i C supp Pi. 

Since Ei C S^+i, using the so-called knot insertion formula, all B-splines of level ^ can 
be written as a linear combination of B-splines of level ^ + 1. More precisely, if Pi G Bi, 
then, 

Pi= ^ c^^^PPi)Pi+i, (3) 

/3^+iSC(/3^) 

where the coefficients Cjs^^PPi) are positive, and C{Pi) C Bi^i is the set of children of Pi. 


Tensor-product B-spline bases Let d > 1. In order to define a tensor-product 
d-variate spline function space on hi := [0,1]*^ C we consider p := (pi,p 2 , ... ,Pd) 
the vector of polynomial degrees with respect to each coordinate direction and n : = 
(ni,n 2 ,... ,nd), where rii > pi + 1. For i = 1,2, ... ,d, let Ep.^n, ■= be a 

Pi-open knot vector, i.e.. 


0 = = 


_ Ai) 

^pi+1 


<i. 


(i) 


< 


'Pi+2 — 


< 5“' 


rii 


^ Sni+l 


f(*) _ 1 

Sni+pi+l 


where the two extreme knots are repeated p* -|- 1 times and any internal knot can be 
repeated at most p* -|-1 times. We denote by 5p_n the tensor-product spline space spanned 
by the B-spline basis i3p,n defined as the tensor-product of the univariate B-spline bases 
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, i3(Sp^ , 1 ^). Let Qp^n be tensor-product mesh consisting of the elements 
Q = Ii X ■■■ X Id, where Jj is an element (closed interval) of the z-th univariate mesh 
X(Hp,,nJ, for z = 

As we did for the univariate case, we consider a given sequence of tensor- 

product d-variate spline spaces such that 

5o C 5i C ^2 C ^3 c ..., (4) 

with the corresponding tensor-product B-spline bases denoted by 

Bo, 6 ^, 62 , 63 ,..., 

respectively. More precisely, if p = {pi,P 2 , ■ ■ ■ ,Pd) is given, for i G Nq) <Si ■= ‘5p,n£ is 
the tensor-product spline space and Bi := Bp ^i is the corresponding B-spline basis, for 
some = {n\,n 2 , ■ ■ ■ ,n^d)- order to guarantee (4), we assume that if ^ is a knot in 
Sp. with multiplicity m, then ^ is also a knot in ^e+i with multiplicity at least m, 
for z = 1,..., d and i G Nq- Furthermore, we denote by Qi := Qp,n^ the tensor-product 
mesh and we say that Q G is a cell of level 1. We state some well-known properties of 
the B-spline basis functions that will be useful in this presentation [dVOl, S07]: 

• Local linear independence. For any nonempty open set O <Z Ll, the functions in Bi 
that do not vanish identically on O, are linearly independent on O. 

• Positive partition of unity. The B-spline basis functions of level I form a partition 
of the unity on i.e., 

'^(5 = 1, on (5) 

y&Bi 

• Two-scale relation between consecutive levels. The B-splines of level i can be written 
as a linear combination of B-splines of level £ -f 1. More precisely, 

Pi = ^ C/3^+i(/5£)/^£+1, y Pi ^ Bi, (6) 

A+iSC(/3f) 

where the coefficients cy^^.^ {Pi) are positive and can be computed using the corre¬ 
sponding coefficients in the univariate two-scale relation (3) and Kronecker prod¬ 
ucts. Here, C{Pi) is the set of children of Pi, and we say that G is a 
child of Pi if the Tth univariate B-spline which dehnes /?£+! is a child of the z-th 
univariate B-spline dehning Pi, for each coordinate direction z = 1,..., d. 

Remark 2.2. Notice that if we dehne cy^^PPi) := 0 when is not a child of Pi, then 
equation (6) can be written as 

Pt= ^ cid^^pPi)Pi+i, ypieBi. (7) 

/3^+ieBf+i 

In particular, we remark that 

C(Pi) = {Pi+i G Bi+i I C/3^+i(/?£) > 0} C {Pi+i G Bi+i I suppPi+i C supp/34- (8) 
Finally, we also consider the set of parents of Pi+i G Bi+i given by 

V{Pi+i) := {Pi G Bi I Pi+i is a child of Pi). 
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3 Hierarchical B-spline basis 


Definition 3.1. If n G N, we say that := {flo, • • •, ^n} is a hierarchy of subdo¬ 
mains of 12 of depth n if 

(i) 12^ is the union of cells of level 1 — 1, for 1 = 1, 2 ,..., n. 

(ii) 12 = 12o D 12i D ■ ■ ■ D 12„_i D 12„ = 0. 


We now define the hierarchical B-spline basis l-L = in the following recursive 

way 


f no := Bo, 

< ni+i := {(3 e nil supp/3 <fL 12f+i} u {/3 e | supp/3 C 12 £+i}, 1 = 

\n-.= nn-x. 

\iVi\=Ui \ He+i, then 

Vi = {13 ene \ supp/3 C 12£+i} = {/3 e Bi\ supp/3 C ll^+i}. 
Notice that in order to get from we replace the set by 

{{3 e I supp/3 C fl£+i}. 


0,..., n — 2. 

(9) 


Moreover, it is easy to check that 

n—1 

G Bi I supp (3 Gfli A supp (3 ^ 12f+i}. (10) 

1=0 

We say that (3 is active if /3 G The corresponding underlying mesh Q = Q(12n) is 
given by 

n—\ 

Q ■= UW e Qi \ Q A Q (f. 12f+i}, (11) 

£=0 

and we say that Q is an active cell is Q E Q, or that Q is an active cell of level i if 
Q E Qn Qi- 

Our definition of the hierarchical B-spline basis n is slightly different from the one 
given in [K98, VGJSll], because they consider the domains 12^, the cells Q and the 
function supports as open sets. That definition does not allow to rebuild the hierarchical 
space uniquely from the only knowledge of the hierarchical mesh (i.e., the active cells of 
each level); see for example [K98, Figure 2.5] or [K97, Figure 2]. 

We know (cf. [VGJSll]) that functions in "H constitutes a linearly independent set 
and that 

So = span So C spanH. 

Unlike the B-spline bases Be for tensor-product spline spaces, the hierarchical B-spline 
basis n does not constitute a partition of the unity. Instead, we can prove the following 
result. 


6 


Lemma 3.2. (Partition of unity in Pi). Let Pi he the hierarchical B-spline basis associated 
to the hierarchy of subdomains of depth n, fin ■= , ^n}- Let := 1 for all 

/So £ ^0 o-nd 

ayt+i ■= X] V/?£+i e S^+i, supp/?£+i C (12) 

he&Bt 

supp /3 £CO£_|_i 


for £ = 0,1,..., n — 2. Then, 


ayfd = 1, on f2. 

yen 

Remark 3.3. In view of the linear independence of functions in Pi, we have that the set 
{ay}y^T-L is uniquely determined. On the other hand, we remark that the dehnition (12) 
depends on the hierarchy of subdomains and that the coefficients ay are dehned not 
only for /S G "H but also for (3 G ^ I supp fii C fLi\. 

Proof. By (5) we have that Yhy^Uo = 1’ 12- Assume now that for a hxed i 

{0 < i <n — 2), Yhy^Uz Thus, using (12) and (7) we have that 


Y1 

y&'Ht+z 



^ (^yi+iZi+i 

y&He 

/3^+iSS^+i 

supp /3 (z!0£_|_i 

suppft+iCO^+i 

Z 

ay (3 + ^ ay,,(3£ 

y&iz 

yt&Bz 

supp/3(z;0^+i 

supp/3^cn^+i 

^ ayl3. 


y&ii 



Thus, E/3eW<+i 

Finally, the proof is complete regarding that Pi = Pin-i- □ 

Notice that in the last lemma ay > 0, for all (3 & Pi. The next result characterizes 
the functions (3 whose weight ay is equal to zero. 

Theorem 3.4. Let G iS^+i he such that supp fti+i C for some i = 0,1,... ,n — 
2. Then, the following statements are equivalent: 

(^) = 0. 

(ii) 'P(3ee V{Ze+i), (ay,, > 0 A supp/S^ C ^ supp/?^ O^+i). 

Proof. From (12), we have that ay^^.^ = 0 if and only if for all /3£ G Bi, with supp /3i C O^+i, 

ay^ = 0 or (/?£)= 0, 

or, equivalently, if and only if, all parents (3£ of /S^+i with ay,, > 0, satisfy supp/d^ Z 

Q 
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Notice that there is no function {3^ G Bq satisfying = 0. Now, the last theorem 
implies that = 0 for (3i G Bi with supp /9i C i if and only if all its parents are active. 
Roughly speaking, ffi is too narrow around supp/3i. Notice that, as soon as one of the 
/Si’s parents is deactivated, will become positive in the new conhguration. 

We conclude this section with the following result, which states that each deactivated 
B-spline of level ^ can be written as a linear combination of functions in the hierarchical 
basis of the subsequent levels, i.e., £ + 1,..., n — 1. 

Lemma 3.5. Let Ti be the hierarchical B-spline basis associated to the hierarchy of sub- 
domains of depth n, f2„ := {flo, ■ ■ ■, f^n}- Then, 

Vi C span I "H n ^ j , (13) 

V k=i+i ) 

for £ = 0,1,..., n — 2. 

Proof. Notice that (13) holds for £ = n — 2 due to (7). Let us assume that (13) holds 
for some £, and prove that it holds for i — 1. Let G P^-i. Since [3i-i G and 
supp/S^-i C Ll£, we have that 

f^e-l = ^ Ci3^{f3£_i)l3i = ^ C/3^ (/?£_! )/?£ + ^ C/3^ (/?£_! )/?£. 

0i&Be ye&Bi /Si&HnBi 

supp/3^cn^ supp^^ci^^+1 

Thus, (3i-i G span (PL fl ljfc=£ ^k), which concludes the proof. □ 


4 Approximation properties and quasi-interpolation 


Let Ptn '■= ^ 1 , • • • 5 ^n} be a hierarchy of subdomains of fl of depth n and let Q be 

the hierarchical mesh given by (11). Let PL be the hierarchical basis dehned in (10). 

A multiscale quasi-interpolant operator has been introduced by Kraft [K98], where 
pointwise approximation estimates were established, for the case fl = where there is 
no boundary. There, the (inhnite) knot vector in each direction considered for building 
the initial tensor-product space Sq was the set of the integer numbers Z, and then dyadic 
rehnement is performed to obtain the subsequent levels. In particular, we extend the 
Kraft construction to the case of open knot vectors with possible multiple internal knots. 
Moreover, we present a multiscale quasi-interpolant operator which provides suitable local 
approximation orders in L'^-norm, for any 1 < g < oo. 

For each coordinate direction i = 1,..., d, we assume that the sequence of knot vectors 
{Spi,nf}teNo satishes: 

• The sequence of meshes is locally quasi-uniform with parameter 6** > 0. 

We let 0 := maxj=i . 


The meshsize in each direction is at least halved when moving from a level to the 
next one, i.e.. 


h 




< -h 


( 14 ) 


for all £ G Ko, where hi^i is the maximum length of the intervals in X(Sp, „^). 
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Remark 4.1. Notice that if the initial knot vectors for each direction, i.e., ,^o, 

..., „o, are chosen arbitrarily and we perform dyadic rehnement for obtaining the 

subsequent levels, we will obtain a sequence of meshes which satishes the two conditions 
just stated. 

In order to dehne the multiscale quasi-interpolant operator in the hierarchical space, 
we need to introduce hrst some local quasi-interpolant operators P^, for £ = 0,1,..., n — 1, 
satisfying certain suitable properties. These last operators can be dehned using the ideas 
from [LLMOl], where each operator is dehned using a underlying local approximation 
method. 


A local approximation method. We recall that p := {pi,p 2 ,... ,Pd) and denote 
by Pp the set of tensor-product polynomials with degree at most Pi in the coordinate 
direction Xi, for i = 1,2,... ,d. Let N := dimPp = -|- 1). 

Let £ G No he arbitrary and hxed. For Q & given, we consider the basis Bq : = 
{Pi ,..., P^} of Pp, consisting of the B-spline basis functions in Bi which are nonzero 
on Q. 

Let LIq : L^{Q) —)■ Pp be the L^-projection operator dehned by 


'Q 


if-UQf)g = 0, V^^ePp. 


(15) 


Notice that 


N 


= V/6n(Q), (16) 

i=l 

where A«(/) := (\? {f is the solution of the linear system 


where 






'Q 


Mq^ = Fq, 

and Fq=( f fP? 




dVxI 


Q . Til 


Since fig preserves polynomials in Pp, we have that {A^ : L^{Q) —)■ M | f = 1,..., N} 
is a dual basis for Bq in the the sense that 


A®(0?) = 


1, if i = j 

0 , 


i,j = 1,...,N. 


(17) 


As a consequence of the L°°-local stability of the B-spline basis we have the following 
result (cf. [BGGS15]). 

Lemma 4.2. Let q he such that 1 < q < oo. Let Q & Qi and let fig : L^{Q) —)■ Pp be 
the L"^-projection operator defined by (15). Then, there exists a constant Cp g > 0 which 
depends only on p and 9 such that 

l|A«(/)||oa < Cp,9|Ori||/||„|^,,, V/ e LHQ), 

where \^{f) = (Af (/),..., A®(/))^ are the coefficients o/ng(/) with respect to the local 
basis Bq as given (16). 
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A locally supported dual basis. For £ = 0,1,..., n — 1, we define oji as the union of 
the elements of level £ whose support extension is contained in i.e., 

'■= U 

where Q denotes the support extension of Q, given by Q = supp /3. In other 

P&Bt 
supp /3 dQ 

words, bJi consists of the elements of level I where the full tensor-product space of level I 
can be exactly represented in the hierarchical space. 

Let 

■.= {13 e Bi\3Q e Qi : Q c. supp {3 n ue}, (18) 

and 

:= spanSf,^,. 

We remark that G {(3 ^ Bi\ supp/? C but in general, Bi^^^ £ £ 

Bi I supp/9 C see Figure 2 below. The goal of this paragraph is to dehne a dual 
basis for the multivariate B-spline basis Be^^Ji, he., a set of linear functionals 

{A/3 : L^(f2) —)■ M I /3 G 

{ 1, % — J 

’ , for all Pi, Pj G Bii^f,. We will use the technique 

0, if i ^ j 

presented in [LLMOl] together with the local L^-projection dehned in (15). Roughly 
speaking, we dehne the functional A/? as a local projection onto some Qp ^ Qe such that 
Qp C supp/d n (Uf. More precisely, for each P G we choose Qp G Qi such that 

Qp G supp P HoJi and let 

:= A®'. 

where zq = io{P, Qp) with 1 < Zq < -A^ is such that = P on Qp. 

As an immediate consequence of (17) and Lemma 4.2 we have the following result. 

Proposition 4.3. Let {A/? : L^(r2) —)■ M |/3 G Bi^i^^} be the set of linear functionals just 
defined above. Then, the following properties hold: 

(i) Local support: IfQp denotes the element in Qi chosen for the definition of Xp, then 
Xp is supported in Qp, i.e., 

^feUiSi), ^ A^(/) = o. 


pi) Dual basis: For Pi, Pj G XpPPj) 


1, ifi = j 

0 , ifi^J 


(Hi) L^-Stability: If P E and f G L'^{Qp), 


\Xp{f)\<C^MQgrh\f\\L.iQ,). 


(19) 
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Localized quasi-interpolant operators in the tensor-product spaces. Now, we 
are in position of defining a quasi-interpolant operator for each level i, using the dual 
bases introduced in the previous paragraph. For i = 0,1,... ,n — 1, let Pi : —)■ 

span C Si be given by 


Pef-= V/eL»(fi). (20) 

The next result summarizes the main properties of Pi. 

Theorem 4.4. For ^ = 0,1,... ,n — 1, let Pi be the operator given by (20). Then, the 
following properties hold: 

(i) Pi preserves splines in i.e., Pis = s, for all s G 
(a) Pi is supported in ui, i.e., 

WfeL^in), ^ Pif^O. (21) 

(Hi) For all s G Si, 

PiS = s, on ooi. 

(iv) Stablility.- The guasi-interpolant operator Pi satisfies 

\\Pif\\Li{nt)^Cs\\f\\L<i{uje), V/GL^(a;£), (22) 

where the constant Cs > 0 only depends on p and 6 . 

(v) Approximation.- Let s := (si, S 2 ,..., Sd) be such that Si <pi + 1, for i = 1,2,... ,d. 

For f G L®(fi) := {g G Tiioc(^) I Dl\g G T'?(fi), 0 < < s^, i = 1,..., d}, 

d 

Wf - ft/lli.K) <CAYi hZWDSJh’ta,), (23) 

i=l 

where the constant Ca > 0 depends on d, s, p and 6 . 

Proof, (i) This is an immediate consequence of Proposition 4.3 (ii). 

(ii) This follows from (19) and (20). 

(hi) This is a consequence of (i) and (ii). 

(iv) Let Q E Qi such that Q C fli. Taking into account (19) and (5) we have that 
\Pif\< max |A^(/)| < C|Qr^||/|| 

supp /3 dQ 


for all / G L'^{ui), where Q 
and 6 . Then, 


= supp /3. Here, the constant C > 0 depends on p 

supp /3 dQ 

l|a/lll.(Q) < C||/||„,gn„,,. (24) 
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Now, (22) follows from the last equation. 

(v) Let Q E Qi such that Q (Z U£. By results on multidimensional Taylor expansions, 
there exists pq G Pp such that 

d 

Wf-PqUm < CT^fc?.l|BS/||„,Q,. (25) 

i=l 

where the constant Ct > 0 only depends on d, s, p and 9. Taking into account (24), (iii) 
and (25) we have that 


11 / - Pe.f\\Li(Q) < 11 / - VoWh^iQ) + WVq - Pe.f\\Li(Q) 

= 11/ -PQlkno) + \\Pti.PQ - f)\\L^{Q) 
— P P')\\f ~ Pq\\li{Q) 

d 

<{l + C)CTY,l^t\K‘J\\LHQy 


Now, (23) follows. □ 

A multiscale quasi-interpolant operator. Let Pe be the operator given by (20), for 
each £ = 0,1,..., n — 1. We dehne If : L^{Q,) —)■ spanTf by 

(Llo ;= Po, 

< n^+i := + P£+i(id-nf), i = 0,...,n-2. (26) 

[n ;= n„_i. 

Remark 4.5. Notice that Theorem 4.4 (i) implies that Pos = s, for all s G Sq. Thus, as 
an immediate consequence of the dehnition of If given in (26), we have that 

Lis = s. Vs G iSo, 

i.e.. If preserves splines in the initial level and in particular, tensor-product polynomials 
in Pp. 

Remark 4.6. Taking into account the dehnitions of If and Pi given by (26) and (20), 
respectively, we have that 


71—1 f n —1 

n ; L'^{Q) — > ;= < s^ | s^ G £ = 0,1,..., n — 1 

r=o I e=o 

When restricting to a hxed ug, the following result allows to write the multiscale 
quasi-interpolant If in terms of the operators Pi, Pt+i,- ■ ■ ,Pn-i- 

Theorem 4.7. If 

Un-l C UJn-2 C ' ' ' C Ci;2 C (Ui C Uq, 



12 


then 


n—1 

= Pif + Pkif - Pk-if), onui, (£ = 0,1,... ,n - 1), 

k=e+i 

for f e L3((]). 

Before proving Theorem 4.7 we state the following elementary result. 

Lemma 4.8. For f G 

= Pif, onuji, (£ = 0,1,... ,n - 1). 

Proof. Let 1 < £ < n — 1 and / G L^{VL). Since G C by Theorem 4.4 

(iii), we have that Pgjl^-if = Ti^^if, on uji. Now, the dehnition of yields II^/ = 
Ili_if + Pi{f - Ili_if) = Pif, on Ui. □ 

Proof of Theorem 4-7. From the dehnition of If given by (26) we have that 

n—1 

n = Pfc(id—nfc_i). 

Now, since Uk C cok-i, using that Pk is supported in ojk (cf. (21)) and Lemma 4.8 we have 
that 

n—1 n—1 

Pfc(id-nfc_i)= Pfc(id-Pfc-i), 

k=e+l k=i+l 

which concludes the proof. □ 

Now, we state and prove the main result of this section. 

Theorem 4.9 (Quasi-interpolation in hierarchical spline spaces). Assume that 

OJn-l c ^n-2 C ■ ■ ■ C a ;2 c CJi C OJq. 

Let s := (si, S 2 , ■ • •, Sd) he such that 1 < Si < pi + 1, for i = 1, 2,..., d. // If : L'^(ff) —)■ 
spanP is the multiscale quasi-interpolant given by (26), then, 

d 

11/ - n/||„,„,, < CVl + 2Cs) YiK-iimthna,). C = 0,1,,,. .n - 1. 

i=l 

for f G Ll{n). 
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Proof. Let / G and let i be such that 0 < £ < n—1. Then, using Theorem 4.7, (22), 

that cjfc C cjfc-i, (23) and (14), we have that 


n-l 


n—1 

= ll/-fi/llMM+ E \\Pdf - Pi,-J)\\Lna,) 

k=e+i 

n—1 

< 11/ - PifllLiiuje) + Cs E 11/ - A.-./iu.(^.) 

k=i^l 

n—1 

< 11/ - n/ll«K) + Cs E 11/ - a-i/llw(„,_.) 


k=£+l 


<CA{j2hZ\\D 


n—1 d 


. i=l 


XiJ II 


n—1 


+cs E E'‘J/i.<ii^s/ii 




A:=£+l i=l 


<CsE i+cs E 


2 = 1 
d 


k=i+l 
oo 


2(A:-1-Qsi 




< Cs E 1 +Cs E ^ La II ®t/ii «(«<) 


2 = 1 
d 




2=1 


/c=0 


Cc 


1 - 2-^» 




□ 

Remark 4.10. If the hierarchy of subdomains Pin '■= {flo) • • • 5 ^n} satishes 

C i = (27) 

we say that the mesh Q is strictly admissible (of class 2, cf. [BuGilS]). In particular, 
if the mesh if strictly admissible, in view of [GJS14, Proposition 20], we have that the 
functions in the truncated basis [GJS12] which take non-zero values on any active cell 
belong to at most two different levels. 

Notice that if a mesh is strictly admissible then satishes 

OJn-l C OJn-2 C • • • C Ci;2 C CCl C OJq. (28) 

On the other hand, in Figure 1 we show some non strictly admissible meshes which 
satify (28). 

We conclude this section by applying Theorem 4.9 to the case of strictly admissible 
meshes. More precisely, we obtain optimal rates of convergence in each level of the 
hierarchical mesh when considering the asymptotic behavior (cf. [SM14, Example 2]). 
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Figure 1: Some examples of four-level meshes for splines of maximum smoothness; we consider 
quadratics {p = 2) on the left and in the middle, and cubics {p = 3) on the right. In all cases, 
the meshes are not strictly admissible but they satisfy (28). The domains wi, ui 2 and wa are 
highlighted in grey from the lightest to the darkest. 


Corollary 4.11. Assume that each level is obtained by dyadic refinement of the elements 
of the previous one (see Remark 4.1). If the mesh is strictly admissible (cf. (27)), then, 

d 

11/ - n/||„(n,, < ft" ||C2/||„,n,_.|, (< = 1,... .n), 

i=l 

for all f G where the constant C > 0 depends only on d, s and p. 


5 A new easier hierarchical spline space 

Assume that we have already computed the set Hi for given £ (cf. (9)). Now, in or¬ 
der to compute Hi+i we need to select the new B-splines to be added, i.e., {(3 G 
Bi+i \ supp/d C Once we know the B-spline basis functions to be deactivated 

Vi = {fi ^ Bi \ supp (3 C Of+i}, notice that it is not enough replacing the functions in Vi 
by their children, because in general, 

U c{fi) £ {/3 G Bi+I I supp/? C 0£+i}. 

y&Vt 

In Figure 2 we show some examples of this situation. This observation suggests a 
simplihed way of selecting B-splines at different levels which consists in adding solely the 
children of the deactivated functions. Doing so, we obtain a new hierarchical space whose 
basis that we call H = Hifiln) is dehned as follows: 

'Ho-.= Bo, 

'Hi+i\= {(3 &iii\svippl3 (fiVLi+i}\J y C(/?), £ = 0,...,n-2. 

supp/JcO^+i 

n:=nn-i. 
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Figure 2: Some examples of two-level meshes for cubic splines {p = 3) of maximum smoothness. 
In both cases, the highlighted B-splines of level 1 have support included in Qi, but they are not 
children of any deactivated B-spline of level 0. 


In this case, iiVi \ Hi+i, we have that 

Vi = {[3 eHi] supp/? c c {/3 e I supp/? c (30) 

but now, we can get T-Le+i from Tic by replacing the B-splines in Vg by their children. 

Thus, it seems that building the basis H is easier than the basis H. In particular, 
there is no need of traversing the mesh in order to identify the B-splines to add in each 
recursive step of (29). However, as an immediate consequence of the following lemma we 
have that 

ncH, 

and therefore, in general, span?/ can be smaller than span?/. 

Lemma 5.1. 

HiCHt, / = 0,l,...,n-l. (31) 

Proof. Notice that (31) holds for / = 0 due to PLq = PLq = Bq. Now, using mathematical 
induction and taking into account (30) and (8) the proof can be completed. □ 

Thus, since that PL C PL, when considering the basis PL instead PL for discretizations in 
isogemetric methods, it will be important to understand which functions we are discarding 
from the basis and the properties of the space span?/. Regarding the set of coefficients 
for the partition of unity in PL (cf. Lemma 3.2) and using Theorem 3.4, we can 
establish the following characterization for functions in PL. 

Theorem 5.2. For i = 0,1,... ,n — 1, 

PLi = {P eHe\ap> 0}. (32) 


In particular, 

PL = {/3 G ?/ I ci/3 > 0}. 
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Proof. Since that ap = 1, for all {3 G Sq, we have that (32) holds for i = 0. Assume now 
that (32) holds for some i and prove that it holds for £ + 1. 

Let (3 G 'Hi+i- If /3 G n using the induction hypothesis we have that 

00 > 0. On the other hand, if /? G PLi+i \ PLi, there exists f3£ E PLi H V{(3) such that 
supp/9^ C and Theorem 3.4 yields O/? > 0. Thus, by Lemma 5.1 we have that 
PL£+i C {13 e PLe+i I a ;3 > 0}. 

Now, let (3 G PLe+i satisfying 0^3 > 0. If /3 G PLe+iPPLi, using the induction hypothesis 
we have that {3 G PLe+i. On the other hand, if /3 E PL^+i \ PLe, we have that (3 G 
and supp/3 C Thus, using Theorem 3.4 we have that there exists (3(, G V{I3) such 

that > 0 and supp/d^ C The induction hypotesis now implies that /?£ G Pi^ and 

therefore, {3 G PLi+i. In consequence, {/3 G PLi+i | > 0} C PLi+i, which concludes the 

proof. □ 

Notice that the functions in PL are linearly independent because PL G PL. On the other 
hand, from (29) it follows that 


(33) 


span C span'Hf+i, £ = 0,1,..., n — 2, 


and therefore, taking into account that PLq = Bo and PLn-i = PL., 


So = spanfdo C spanH. 


Remark 5.3. Since spanido C spandf, we have that tensor-product polynomials in Pp 
belong to spanTi. 

In the Section 5.1 we study the local approximation properties of the space spandf, 
through multiscale quasi-interpolant operators. 

Finally, taking into account Theorems 5.2 and 3.4 we can prove the analogous result 
of Lemma 3.5 when considering the basis PL. 

Lemma 5.4. Let PL he the hierarchical B-spline basis defined by (29) associated to the 
hierarchy of subdomains of depth n, fin ■= {flo, • • •, f^n}- Then, 



(34) 


for £ = 0,1,..., n — 2. 

Proof. Notice that (34) holds for £ = u — 2 due to (6). Let us assume that (34) holds 
for some i, and prove that it holds for i — 1. Let G Vi-i. Since G and 
supp/d£_i C fl£, we have that 


fie-l — ^2 Cyfil3£-i)l3£ — ^ C/3^(/?£_!)/?£-1- ^ Cyfil3£-i)l3£. 


/3^£C(/3£_i) 




0l£C{0(-i)\'Di 


Thus, f3£-i G span (pLH Ufc=£ ), which concludes the proof. 


□ 
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5.1 Quasi-interpolation and local approximation properties in 

span?/ 

In this section we assume that each level is obtained by dyadic rehnement of the elements 
of the previous one (see Remark 4.1). The following related auxiliary technical result will 
be useful. The proof is presented in the Appendix 1. 

Lemma 5.5. Let be the set of B-splines defined in (18), for I = 0,1,... ,n — 1. 
Then, 

C IJ £ = 0,l,...,n-2. 

P&Bt 

supp 

This lemma allows us to prove the following proposition, which, together with the 
results presented in the previous section, show that the space span Li is rich enough, and 
in particular contains all the local spaces 

Proposition 5.6. Assume that 

OJn-l c OJn-2 C • • • C Ci ;2 C (Ui C OJq. (35) 

Then, 

£ = 0,1,... ,n - 1. 

Proof. Taking into account the definition of PLi given in (9) and the characterization of 
idi in (32), we have that {fi & Bi \ supp /? C A > 0} C He, for £ = 0,1,..., n — 1. 
On the other hand, for fi G Be^^^, we have that supp fi (Z fie and thus, ap is well-dehned 
(cf. (12)). Thus, it will be enough to prove that 

Be,ui Z {(i eBe\ag> B], £ = 0,1,..., n - 1. (36) 

Notice that (36) holds for £ = 0 due to = 1 > 0, for all fio E Bq = Bq^ujq- Now, 

assume that (36) holds for some f. Let fie+i ^ Be+i^ut+i- In view of Lemma 5.5, there 
exists fie G 'B{fie+i) such that supp/3^ C O^+i. Moreover, taking into account (35), the 
dehnition of Be+i^u^+i implies that there exists Qe+i G Qe+i such that 

Qe+i Z supp fie+i O ue+i C supp fie O oje, 

which in turn yields fie G Be,Lo^. Finally, taking into account the induction hypotesis, we 
have that > 0 and now using Theorem 3.4 we conclude that >0. □ 

In view of Remark 4.6, the immediate consequence of Proposition 5.6 is that the 
multiscale quasi-interpolant operator dehned in Section 4 does construct interpolating 
functions belonging to PL and not only to PL. Thus, the approximation estimates from 
Theorem 4.9 apply verbatim to PL. We express this fact in the following simple corollary. 

Corollary 5.7. Assume that 

tOn-l C UJn-2 C • • • C Ci ;2 C (Ui C LUq. 

Let 1 < q < oo and 11 : L^{fl) —)■ spanH he the multiscale quasi-interpolant operator 
defined in (26). Then, 

n : Lfifl) -E spanPL. 
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5.2 Refinement of hierarchical spline spaces 

When thinking of hierarchical splines within a rehnement and an adaptation process, 
it is very important to have a precise link between the enlargement of the hierarchy of 
snbdomains ■ ■ ■; and the rehnement of the corresponding hierarchical 

space. This issne has been addressed for classical hierarchical splines in [GJS14] and here 
we address it for T-L. 

Definition 5.8. Let fin ■= (ffo) ! ^n} and := (ffo) • ■ ■) ^n+i} b® 

hierarchies of snbdomains of hi of depth (at most) n and n + 1, respectively. We say that 
is an enlargement of fin if 

fliCfl}, i=l,2,...,n. 

Let be an enlargement of fin- Now, the corresponding hierarchical B-spline 

basis and refined mesh Q* are given by 

n 

T-L* := l^{/d e Bi I supp fi (Z fl*g A snpp fi (fi 
£=0 

and 

n 

Q* ■= Uw eQe \ Qzfil A 
£=0 

In [GJS14] has been proved that 

spanH C spanH*. (37) 

Let denote the sequence of coefficients (with respect to the hierarchy Li)i_|_i) 

given by Lemma 3.2. Thus, we have that 

^ a*^l3 = 1, on n, 
g&H* 

and thus, we can consider 

W := {fieW \a}> 0}. 

The following theorem establishes the analogous of (37) when considering the basis H. 

Theorem 5.9. 

span 77 C span 77*. 

In order to prove this result we need the following auxiliary lemma. 

Lemma 5.10. If (3 E and supp fi C fli, for i = 0,1,... ,n — 1, then 

®/3 ^ cbg- 
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Proof. The assertion holds for £ = 0 due to = 1 for all {3 G Bq. Now, assume 

that for some ^ we have that 

> Oft, for Pi G Bi, such that supp/3^ C fli. (38) 

Let Pi+i G Bi+i such that supp/S^+i C Since C using (38) and the 

dehnitions of and we have that 

^h+i= Y1 > Y1 aAC;3,+i(/5^) = a/3,+1 • 

/3,eB, ft SB, ft SB, 

suppftcO|+i suppftcO,+i suppftcn,+i 


□ 

Proof of Theorem 5.9. Let P E Pi and let £ be such that P G Bi. Since supp P <zPLi C. 
we have that P G PL}. On the other hand, Lemma 5.10 implies that a*^> > 0, and (32) 

yields P E PL}. Finally, taking into account (33) we have that P E span?i*. □ 

6 Concluding remarks 

In this paper, after studying the approximation properties of hierarchical splines as de- 
hned in [K98], we propose a new hierarchical spline space, through a construction of a 
set of basis functions named PL, that enjoys several properties and may be considered as 
a valuable alternative to truncated hierarchical splines [GJS12]. We can summarize and 
comment upon our results as follows: 

• The basis that we construct simplihes the implementation and data structures 
needed to carry hierarchical splines because the rehnement can be performed through 
the parent-children relations between B-splines. Unlike the classical hierarchical 
space, where algorithms traversing the mesh are needed to identify the new active 
B-splines, we just add children of already active B-splines. Moreover, we believe 
that our construction can be suitably used in conjunction with function-based er¬ 
ror indicators, i.e., error indicators that mark functions (and not elements) to be 
rehned: in our framework, a marked function would be simply replaced by some of 
its children. These aspects are studied in the forthcoming paper [BuGal5]. 

• The weighted basis {a^/? | P G PL} constitutes a convex partition of unity and has 
the advantage of preserving simple basis function supports (only hypercubes) and 
also, in principle, ask for the use of simple spline evaluation formulae. 

• We have extended the multiscale quasi-interpolant proposed in [K98] to the cases of 
general and open knot vectors and we have provided local approximation estimates 
in L'^-norms, for 1 < g < cx). This interpolant is built in a multiscale fashion, it 
is not a projector in general and is not based on dual functionals, unlike the one 
presented in [SM14] based on THB-splines. On the other hand, it verihes optimal 
local approximation estimates. 
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Appendix 1 

Now we present the proof of Lemma 5.5. Here, if A and B are two knot vectors (i.e., 
sequences), we say that A C B if A is a subsequence of B. On the other hand, given two 
arbitrary sequences A and B, we denote hy AA B the largest subsequence of A and B. 

Proof of Lemma 5.5. For clarity of presentation, we consider hrst the univariate case 
d = 1. Let p = Pi be the polynomial degree and let ^ be hxed satisfying 0 < £ < n — 2. 
Let and be the open knot vectors associated to the spline spaces Sg, and 
respectively. Let 

Let /df +1 G Then, there exists Q G Qi+i such that Q C supp /?£+! and Q C cj^+i. 

Thus, there exists k such that Q = [^^,^^+ 1 ]- By the definition of we have that 

Q [^k—pt ^k+p+l\ P 

Notice that 

■ {^k—py ■ ■ ■ )C^A:+p+l} 

consists of 2p + 2 consecutive knots in S^+i. Since is obtained from Ei by dyadic 
rehnement, we have that #(Sq fl > p + 1. If O = p + 1, it is easy to 

check that we can add one knot G such that #((Sg U {.^}) fl = p + 2 and 
{.^ G (Sg U {.^}) n C H^+i. Therefore, there exists Sk c such that 

#s^g > P + 2, n s^) c (Sg n s,) c s^g, {f g sy c (39) 

Let r := fl S^) and notice that r < p + 1. We consider two cases: 

(i) minS; 3 ^^^ or maxS^^_^j matches a knot in By (39), there exists (di G Be, such 
that n Sf C C Sk. Thus, (3e+i e C{(de) and supp/3^ C H^+i. 

(ii) Neither minS^^_^j nor maxS;?^^^ match knots in Since = p + 2, in this 

case r < p. Thus, there exists (de G Be such that fl C C S^. Again, 
jde+i e C{lde) and supp /de C H^+i. 

Finally, for the multivariate case d > 1, we can apply this argument in each coordinate 
direction. □ 
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